Statistical ensembles in Hamiltonian formulation of hybrid quantum-classical systems 
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General statistical ensembles in the Hamiltonian formulation of hybrid quantum-classical systems 
are analyzed. It is argued that arbitrary probability densities on the hybrid phase space must be 
considered as the class of possible physically distinguishable statistical ensembles of hybrid systems. 
Nevertheless, statistical operators associated with the hybrid system and with the quantum sub- 
system can be consistently defined. Dynamical equations for the statistical operators representing 
the mixed states of the hybrid system and its quantum subsystem are derived and analyzed. In 
particular, these equations irreducibly depend on the total probability density on the hybrid phase 
space. 
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Hybrid quantum-classical (QC) systems are neither 
quantum nor classical. There is no unique generally ac- 
cepted theory of the QC systems, primarily because each 
of the suggested theories has some unexpected or con- 
troversial features [H-Q- Current technologies are suffi- 
ciently developed to enable experimental studies of the 
interaction between typically quantum and typically clas- 
sical objects @. The choice of the proper hybrid theory 
will ultimately depend on the experimental tests, but 
such experiments require detailed preliminary theoreti- 
cal models. It is plausible to expect that the interac- 
tion between the classical and the quantum subsystems 
might introduce features that are not present neither in 
the quantum nor in the classical subsystems without in- 
teraction (see for example [l| ) • In this communication we 
shall explore some features of the hybrid system theory 
which is formulated usingthe framework of Hamiltonian 
dynamical systems [l], [1,0] ■ 

The Hamiltonian hybrid theory, as formulated for ex- 
ample in [![, has many of the properties commonly ex- 
pected of a good hybrid theory. However, it also has 
some controversial features concerning the class of math- 
ematical objects that should be interpreted as physical 
variables of the QC system. Also, what should be mathe- 
matical objects that represent the most general statistical 
ensembles of QC systems and the corresponding mixtures 
of the quantum subsystems is not trivially obvious and 
needs a careful discussion. Formally, the problem is in- 
herited from the surplus structure present in the Hamil- 
tonian formulation of the quantum mechanics. It can 
be argued that what must be considered as nonphysical 
in the Hamiltonian formulation of the standard linear 
quantum mechanics might acquire physical meaning for 
the QC system because of the presence of the classical 
subsystems and the quantum-classical interaction. 

In what follows we shall first briefly recapitulate the 
Hamiltonian formulation of quantum mechanics and of 
the hybrid quantum-classical systems. Then we shall dis- 
cuss possible representation of general statistical ensem- 



bles of QC systems within the Hamiltonian formulation. 
Statistical operator for the QC system, as well as con- 
ditional and unconditional mixed states of the quantum 
subsystem corresponding to the general ensembles of the 
QC systems, will be defined and their evolution will be 
discussed. 

Hamiltonian formulation of quantum mechanics 
Schrodinger dynamical equation on a Hilbert space If. 
generates a Hamiltonian dynamical system on an appro- 
priate symplectic manifold @, Q . The real manifold M. , 
associated with the Hilbert space ri in fact has Rieman- 
nian and symplectic structure, provided by the real and 
the imaginary parts of the scalar product, and can be 
viewed as a phase space of a Hamiltonian dynamical sys- 
tem, additionally equipped with the Ricmannian metric 
which reflects its quantum origin. A vector from rl, 
associated with a pure quantum state, is represented by 
the corresponding point in the phase space M. denoted 
by X^ or simply by X. 

Real coordinates {(xj,yj), j = 1,2,...} of a point tp e 
ri = At arc introduced using expansion coefficients {cj, 
j = 1,2,...} in some basis j = 1,2,...} of ri as 
follows 



|V) = Xj = V2Rc Cj, tjj = V2Imcj 



(1) 



The coordinates (xj,yj) represent canonical coordinates 
of a Hamiltonian dynamical system on M. . Consequently, 
the Poisson bracket between two functions F\ and F% on 
M. in the canonical coordinates (xj , yj ) is given by 

r p p 1 _l^ r (dF 1 dF 2 dF 2 dFA 

A quantum observable H is represented by the corre- 
sponding function of the form 



(3) 
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Hamiltonian flows with the Hamilton's function of the 
form (3) generate isometries of the Riemannian metric. 
More general Hamiltonian flows on M. , corresponding to 
the Hamilton's function which are not of the form (3), 
do not generate isometries and do not have the physical 



2 



interpretation of quantum obscrvablcs. It can be seen 
easily that 



{H X ,H 2 } M = -{[H U H 2 ]). 
in 

The Schrodinger evolution equation 

ih\i>) = H\i>) 



(4) 



(5) 



is equivalent to the Hamilton equations on Ai assuming 
the standard form in the canonical coordinates (xj,yj) 

x o = 7Tr> Vi = ~ — i ( 6 ) 



dx-j 



with H given by (3). 

The Hamiltonian formulation of the Schrodinger equa- 
tion on T-L automatically preserves the constraints im- 
posed by the physical equivalence of Hilbert space vec- 
tors. In fact, Hamiltonian formulations based on % and 
on the projective Hilbert space as the space of physical 
states are equivalent. We use the formulation in which 
points of the quantum phase space are identified with the 
vectors from T-L since it is sufficient for our main purpose. 

Mixed states of a quantum system in the Hamiltonian 
formulation 

A quantum state is in general represented by the cor- 
responding density operator p on %. On the other hand 
every positive function p{x 1 y) with unit integral on Ai 
represents a density of some probability theory on Ai. 
Expectation of a function F(x, y) with respect to p(x, y) 
is given by 



F 



P {x,y)F{x,y)dM, 



(7) 



M 



where dM is the Lebesque measure on Ai . The densities 
satisfy Liouville equation on Ai 
d 

g^p(x,y,t) = {H(x,y),p(x,y;t)} M . (8) 

Quantum mechanical average of the observable F in the 
state p, Tr(p-F), is reproduced with the formula (7) using 
any of the probability densities p(x,y;t) with the same 
first moment fixed by the requirement 



P(to) 



p(x,y;t )fL(x,y)dM, 



(9) 



M 



where Tl(x,y) = \ipx,y) (ipx,y\ and the state \4>x,y) € H 
corresponds to the coordinates (x,y) of Ai. Liouville 
evolution of the densities p(x, y; t) yielding the same p(to) 
generates the same von Neumann evolution pit). The 
fact that the quantum mixed state p determines only an 
equivalence class of densities, those with the appropriate 
first moment, is equivalent to the non-uniqueness of the 
expansion of the quantum mixed state in terms of convex 
combinations of pure state projectors. 

Hamiltonian theory of hybrid systems 

Hamiltonian theory of hybrid quantum-classical sys- 
tems can be developed starting from the Hamiltonian 
formulation of a composite quantum system and impos- 
ing a constraint that one of the components is behaving 



as a classical system Q. The result in the macro- limit 
imposed on the classical subsystem turns out to be equiv- 
alent to a Cartesian product of two Hamiltonian systems 
as in One of these Hamiltonian systems corresponds 
to the quantum and one to the classical subsystems of 
the hybrid. However, the interaction between the two 
subsystems has crucial influence on their properties. 

The phase space of the hybrid system Ai is considered 
as a Cartesian product Ai = Ai c x Ai q of the classical 
subsystem phase space Ai c and of the quantum subsys- 
tem phase space Ai q . Denoting the local coordinates 
on the product as {p, q,x,y}, where (jp, q) € Ai c and 
(x, y) <E Ai q one can write the evolution equations of the 
QC system as Hamiltonian dynamical equations on the 
phase space Ai with the Hamilton's function comprised 
of three terms 

H t (p,q,x,y) = H c (p,q)+H q (x,y) + V int (p,q,x,y), (10) 

where H c is the Hamilton's function of the classical sub- 
system, H q (x,y) of the form (3) is the Hamilton's func- 
tion of the quantum subsystem and Vi nt (p,q,x,y) = 

(^x, y \Vint(p,q)\tpx,y), where Vint(p,q) is an operator in 
the Hilbert space of the quantum subsystem which de- 
pends on the classical coordinates (p, q) and describes 
the interaction between the subsystems. The Poisson 
bracket on Ai of arbitrary functions of the local coordi- 
nates (p, q, x, y) is defined as 

'dhdh dhdh 
2i dpi dqi dpi 



if f\ ST ohofr Of 2 a fx 



n^f\dxjdyj dxjdyj 



(11) 



Thus, the Hamiltonian form of the QC dynamics on Ai 
as the phase space reeds 



q = {q,H t }M, 
x = {x, H t }M, 



V = {p,H t } M , 

i) = {y,H t } M , 



(12) 
(13) 



where the Hamilton's function Ht(p,q,x,y) in local co- 
ordinates on Ai is given by (10). 

Statistical ensembles of QC systems and quantum sub- 
systems 

Consider a general probability density p{p, q, x, y) on 
the total hybrid phase space Ai = Ai c x Ai q . There 
is no reason to require such probability density to rep- 
resent a physical quantity or an observable of the QC 
system. If such p{p 1 q, x, y) is a quadratic function of x, y 
then it is equal to the quantum expectation in the cor- 
responding state \ip x ,y) of an operator function f p (p,q), 
i.e., p(p,q,x,y) = (ip x ,y\fp(p,q)\ipx,y), where for each 
fixed p, q the operator f p (p,q) is a statistical operator 
on the Hilbert space of the quantum subsystem. How- 
ever, in general a probability density of an arbitrary form 
p(p, q, x, y) describes a perfectly legitimate statistical en- 
semble of QC systems. In general, following the Hamilto- 
nian formulation of the QC system dynamics, the evolu- 
tion of p(p, q, x, y; t) considered as a statistical ensemble 
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on M. is given by the Liouville equation with the Hamil- 
ton's function (10) and the Poisson bracket (11) 

d 

Q^p(p,q,x,y;t) = {H t (p,q,x,y),p(p,q,x,y;t)} M , (14) 

i.e., p(p(t),q(t),x(t),y(t);t)=const when (p(t),q(t)) and 
(x(t), y(t)) are determined from the Hamilton equations 
(12) and (13), respectively. However, Liouville evolu- 
tion of an ensemble which is at t = to of the form 
(ipx,y\fp(p, q\ to)\ipx,y) will in general result in some prob- 
ability density p(p, q, x, y; t) which is not quadratic in 
x, y, i.e., can not be expressed as expectation of an oper- 
ator. Therefore, it can be argued that the most general 
statistical ensembles of QC systems need to be repre- 
sented by general probability densities p(p, q, x, y; t). We 
shall therefore assume, in accordance with the Hamil- 
tonian theory, that an arbitrary probability density 
p(p, q, x, y; t) describes a statistical ensemble of QC sys- 
tems, that is a mixed state of the hybrid, and that the 
evolution of such mixed states is given by the correspond- 
ing Liouville equation. 

Suppose a QC system is in a general mixed state 
p(j>,q,x,y;t). The density p(p,q,x,y;t) generates a 
unique positive operator valued function (POVF): 



p( t ) = p(P,q,x,y;t)U(x,y)dM. 



(17) 



M 



pip,q,x,y;t)tl(x,y)dM q 



(15) 



which can be called the hybrid statistical operator. 
p(p, q; t) contains less information about the hybrid sys- 
tem state then the density p(p, q, x, y; t), and plays a sec- 
ondary role in the hybrid theory presented here. The cor- 
responding mixed state of the quantum subsystem condi- 
tional on the classical subsystem being in the state (p, q) 
is uniquely represented by 

/3 P ,g(i) = p(p,q;t)/ p(p,q,x,y;t)dM q . (16) 

The unconditional mixed state of the quantum subsystem 
of the hybrid in the state p(p, q, x, y; t) is also uniquely 
obtained as 



At time t the previous formula defines positive, trace 
one operator, i.e., a statistical operator representing the 
mixed state of the quantum subsystem. At any t the 
statistical operator p(t) depends on the value of the sub- 
integral expression at the same time t. The Liouville 
evolution of p(p, q, x, y; t) is certainly continuous in t and 
the dependence on t of p(t) given by (17) is also contin- 
uous. Thus, the formula (17) defines a continuous one- 
parameter family of statistical operators on H.. 

Analogously to the relation (9) valid for a purely quan- 
tum system, many hybrid ensembles represented by dif- 
ferent p(p, q, x, y; t) have the quantum subsystem in the 
same conditional or unconditional mixed state. The cru- 
cial difference between the purely quantum and the hy- 
brid systems is that we have assumed that each different 
pijp, q, x, y; t) describes physically different ensembles of 
QC systems with the quantum subsystem in the same 
mixed state. This will be reflected in the evolution of 
(15) or (17). Recall that all different p(x,y;t) in (9) 
with the same first moment correspond to the physically 
equivalent quantum mixture p(to)-, and generate unique 
von Neumann evolution of p(t) which is obtained from 
the Liouville evolution of any such p(x,y;t). Therefore, 
all such p(x, y\ t) are equivalent in the purely quantum 
case. In the hybrid case, different p(p, q, x, y; t) which 
give the same p(p, q; t ) (or p(i )), as we shall see, gener- 
ate different evolution of p(p, q; t) (or p{t)) and thus must 
be considered as physically different. 

The evolution equation satisfied by p(p, q; t) can 
be obtained from (14) and (15) using partial inte- 
gration over [xj,yj) and the identities dF/dxj = 
((ip*,y\F\j) + (j\P\1>x,v))/y/2, dF/dyj = i((^, y \F\j) - 
(j\F\ip x ^ y ))/y/2, where F = (ip x ,y\F\ip x , y ) . The resulting 
equation is 



J 



dp(p,q;t) = _j_ 

dt ih 



[H q + Vim {p, q),p(jp, q; t)] + {H c {p, q), p(p, q; t)} p q + \ {V int (p, q, x, y), p(p, q, x, y; t)} v g TL(x, y)dM q 



1 p,i 



(18) 



The solution of (18) remains a well defined statistical The equation for the statistical operator of the quantum 
operator on H for all t, which is a desirable property subsystem p(t) follows after the integration over (p, q) 
not shared by some other hybrid system theories [4j, l5j . 



dp(t) 
dt 



^[H q ,p(t)}+^Z [ [V mt (p,q),p(p,q;t)]dM c + [ {H c (p,q),p(p,q;t)} dM c 



M 



{Vintip, q, x, y),p(p, q, x, y; t)} p q IL(x, y)dM. 



(19) 
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The first term on the right side of (19) generate the 
unitary part of the evolution. The last three terms do 
not preserve the norm of p, and are responsible for non- 
unitary effects. Notice that the evolution of p(p, q; t) 
(p(t)) cannot be expressed only in terms of p(p,q;t) 
{p(t)), but irrcducibly involves the probability density 

Observe that, taking different p'(p,q,x,y]t) yielding 
the same p'(p,q;to) = p(p,q;to) via (15) will in general 
generate different p'(p,q;t) ^ p(p,q;t). In other words, 
the states of the quantum subsystem of a hybrid in differ- 
ent states p(p,q,x,y;t) and p'(p,q,x,y;t) might be the 
same at some moment t$, but will inevitably evolve differ- 
ently. This is natural since, the corresponding evolution 
equation (18) for p(p,q;t) must depend on the evolution 
of the entire QC system. In particular, one might ad- 
just the total initial ensemble of the hybrid so that the 
evolution of the quantum subsystem from a fixed initial 
mixture p(to) has different properties, without altering 
the Hamiltonian. Experimental observation of different 
evolutions of the same initial quantum state p(p,q;to), 
obtained from multiple p(p, q, x, y; t) that are different 
functions of {x,y), would provide a confirmation of our 
main assumption concerning the class of physically dis- 
tinguishable ensembles of hybrid QC systems. 

Summary and discussion 

In summary, we have explored some of the conse- 
quences of the assumption that in the Hamiltonian for- 
mulation the set of ensembles of hybrid quantum-classical 
systems is mathematically represented by the space of 
probability densities on the hybrid system phase space. 
Each such ensemble uniquely determines a conditional 
and an unconditional mixed state of the quantum sub- 
system represented by the corresponding density opera- 
tors on the quantum subsystem Hilbert space. Evolution 
of the quantum subsystem mixtures is defined using the 
evolution of the hybrid probability density. Different hy- 
brid ensembles might give the same quantum subsystem 
mixture at some time to, but that quantum mixture ob- 
tained from different hybrid ensembles will evolve differ- 
ently. The evolution equations for the hybrid statistical 
operator p(p, q; t) (18) and for the quantum subsystem 
unconditional mixture p(t) (19) have been derived and 



inevitably involve the full density p(p, q, x, y; t). 

The initial assumption about the mathematical objects 
needed to represent all physically possible ensembles of 
hybrid systems in fact assumes that the ensembles of pos- 
sibly interacting quantum-classical systems are more gen- 
eral then ensembles of quantum-classical systems without 
the interaction between the subsystems. In the later sit- 
uation a set of densities on the hybrid phase space that 
is invariant under the evolution can be chosen to con- 
tain only function that are necessary of quadratic depen- 
dence on the quantum degrees of freedom. However, if 
the quantum and the classical subsystem interact than 
the invariant set of densities is in general the full set 
of probability densities on the hybrid phase space. We 
have shown how such general set of hybrid ensembles 
generates consistently defined quantum mechanical mix- 
tures of the quantum subsystem with the correspond- 
ing evolution equations reflecting the quantum-classical 
interaction. The interaction implies that the evolution 
of the quantum subsystem statistical operator shows ex- 
plicit dependence on the equivalent representations of the 
initial density operator. 

Added note Formula (19) can be considerably sim- 
plified since the last two terms containing integrals over 
classical variables of the Poison brackets are in fact equal 
to zero. However, this does not effect the main conclu- 
sions that the evolution of p(t) depends on the full hy- 
brid density and is non-unitary, since these conclusions 
are consequences of the term with the commutator of the 
interaction V(p,q) and p(p,q). Only if there is no inter- 
action between the classical and the quantum degrees of 
freedom the evolution of the quantum density p(t) is uni- 
tary and independent of its initial expression in terms of 
the full hybrid density. 
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